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ABSTRACT 

We analyse the properties of Mg n absorption systems detected along the sightlines toward GRBs using a sample of 10 GRB afterglow 
spectra obtained with VLT-UVES over the past six years. The signal-to-noise ratio is sufficiently high that we can extend previous 
studies to smaller equivalent widths (typically W T > 0.3 A). Over a pathlength of Az ~ 14, we detect 9 intervening Mg n systems 
with Wy > 1 A and 9 weaker Mgll systems (0.3 < W r < 1.0 A) when about 4 and 7, respectively, are expected from observations of 
QSO sightlines. The number of weak absorbers is similar along GRB and QSO lines of sight, while the number of strong systems is 
larger along GRB lines of sight with a 2cr significance. Using intermediate and low resolution observations reported in the literature, 
we increase the absorption length for strong systems to Az = 31.5 (about twice the path length of previous studies) and find that the 
number density of strong Mgn systems is a factor of 2.1+0.6 higher (about 3tr significance) toward GRBs as compared to QSOs, 
about twice smaller however than previously reported. We divide the sample in three redshift bins and we find that the number density 
of strong Mg n is larger in the low redshift bins. We investigate in detail the properties of strong Mg n systems observed with UVES, 
deriving an estimate of both the H i column density and the associated extinction. Both the estimated dust extinction in strong GRB 
Mg ii systems and the equivalent width distribution are consistent with what is observed for standard QSO systems. We find also that 
the number density of (sub)-DLAs per unit redshift in the UVES sample is probably twice larger than what is expected from QSO 
sightlines which confirms the peculiarity of GRB lines of sight. These results indicate that neither a dust extinction bias nor different 
beam sizes of the sources are viable explanations for the excess. It is still possible that the current sample of GRB lines of sight is 
biased by a subtle gravitational lensing effect. More data and larger samples are needed to test this hypothesis. 
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1. Introduction 

Thanks to their exceptional brightness, and although fading very 
rapidly, Gamma-Ray Burst (GRB) afterglows can be used as 
powerful extragalactic background sources. Since GRBs can be 
detected up to very high redshifts ( |Greiner et al . 2008 1 Kawai 
|et al.|[2006 |Haislip et al. 2006 1 their afterglow spectra can be 
used to study the properties and evolution of galaxies and the 
IGM, similarly to what is traditionally done using QSO spectra. 

Even if the number of available GRB lines of sight (los) is 
much smaller than those of QSOs, it is interesting to compare the 
two types of lines of sight. In particular, |Prochter et al. ( 2006b ) 
found that the number density of strong (rest equivalent width 
W r > 1 A) intervening Mg n absorbers is more than 4 times 
larger along GRB los than what is expected for QSOs over the 
same path length. This result has been derived from a sample of 
14 GRB los and a redshift path of Az = 15.5, and has been con- 
firmed by Sudilovsky et al. ( 2007) 1. Dust extinction bias for QSO 
los, gravitational lensing, contamination from high-velocity sys- 
tems local to the GRB and difference of beam sizes are among 
the possible causes of this discrepancy. All these effects can con- 
tribute to the observed excess, but no convincing explanation has 



( Sudil ovsky et al.|2007 



been found to date for the amplitude of the excess (Prochter et al. 
|2006b[ IFrank et al.||20071 |Porciani et al. 12007) . Similar studies 
'ejos et al.|2007) > have been performed 



on the number of C rv intervening systems. Their results are in 
agreement with QSO statistics. 

Clearly, further investigation of this excess is required. Since 
the reports by Prochter et al. (2006b) (based on a inhomoge- 
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neous mix of spectra from the literature) and the confirmation 
by Sudilovsky et al. ( 2007 1 (based on a homogeneous but lim- 
ited sample of UVES los), several new los have been observed. 
As of June 2008, the number of los with good signal-to-noise 
ratio observed by UVES has increased to 10. We use this sample 
to investigate the excess of strong Mg n absorbers and, thanks to 
the high quality of the data, we also extend the search of systems 
to lower equivalent widths and derive physical properties of the 
absorbing systems. In addition, to increase the redshift path over 
which strong Mg u systems are observed, we consider also a sec- 
ond sample that includes in addition other observations gathered 
from the literature. 

We describe the data in Section|2j identify Mg n systems and 
determine their number density in Section [3] We derive some 
characteristics of strong Mg n systems in Section|4j we estimate 
their HI content in Section|5]and we study peculiar (sub-)DLAs 
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systems detected along the lines of sight in Section[6] We sum- 
marize and conclude in Section [7] 



2. Data 

Our first sample (herafter the UVES sample) includes ten 
GRB afterglows with available follow-up VLT/UVESQ high- 
resolution optical spectroscopy as of June 2008: GRB 021004, 
GRB 050730, GRB050820A, GRB050922C, GRB 060418, 
GRB 060607 A, GRB 071031, GRB 080310, GRB 0803 19B and 
GRB 0804 13 A. All GRBs were detected by the Swift satellite 
( |Gehrels et al.|2004| l, with the exception of GRB 021004, which 
was detected by the High-Energy Transient Explorer (HETE-2) 
satellite ( |Ricker et al.|2003] >. 

UVES observations began on each GRB afterglow with the 
minimum possible time delay. Depending on whether the GRB 
location was immediately observable from Paranal, and whether 
UVES was observing at the time of the GRB explosion, the af- 
terglows were observed in either Rapid-Response Mode (RRM) 
or as Target-of-Opportunity (ToO). A log of the observations is 
given in Table 1 . 



Table 1. GRB sample and log of UVES observations 



GRB 


utF] 


«£| 




ESO Program 


PI 


(yymmdd) 


Swift 


(hn:mm) 


(h) 


ID 




021004 


12:06:13 


13:31 


2.0 


070.A-0599^1 


Fiore 


050730 


19:58:23 


04:09 


1.7 


075.A-060r 


Fiore 


050820A 


06:34:53 


00:33 


1.7 


075.A-0385 


Vreeswijk 


050922C 


19:55:50 


03:33 


1.7 


075.A-0603 


Fiore 


060418 


03:06:08 


00:10 


2.6 


077.D-0661 


Vreeswijk 


060607A 


05:12:13 


00:08 


3.3 


077.D-0661 


Vreeswijk 


071031 


01:06:36 


00:09 


2.6 


080.D-0526 


Vreeswijk 


080310 


08:37:58 


00:13 


1.3 


080.D-0526 


Vreeswijk 


0803 19B 


06:12:49 


00:09 


2.1 


080.D-052#H 


Vreeswijk 


080413A 


02:54:19 


03:42 


2.3 


081.A-085rT 


Vreeswijk 



" UT of trigger by the BAT instrument on-board Swift. Exception: 
GRB 021004, detected by WXM on-board HETE-2. 

b Time delay between the satellite trigger and the start of the first 
UVES exposure: normally a series of spectra is taken. 

Total UVES exposure time including all instrument setups. 

d Also 070.D-0523 (PI: van den Heuvel). 

e Also 080.A-0398 (PI: Fiore). 



The observations were performed with a 1.0" wide slit 
and 2x2 binning, providing a resolving power of R * 48 000 
(FWHM ~ 7 km s 1 ) for a «2 km s 1 pixel sizd^ The UVES 
data were reduced with a customized version of the MIDAS re- 
duction pipeline (Balles ter et al.|20 06). The individual scientific 
exposures were co-added by weighting them according to the in- 
verse of the variance as a function of wavelength and rebinned 
in the heliocentric rest frame. 

Although the UVES sample has a smaller number of los 
compared to the sample used by |Prochter et al. (2006b]) (10 in- 
stead of 14 los), the redshift path of the two samples is simi- 
lar (Az = 13.9 and 15.5 for the UVES and Prochter's samples 
respectively) because of the larger wavelength coverage of the 
UVES spectra. The UVES sample has 6 new los (050922C, 



Dekker et al. 



(20001 



1 UVES is described in ] 

2 Though the minimum guaranteed resolving power of UVES in this 
mode is 43 000, we find that in some cases a higher resolution, up to 
~50000, is achieved in practice, due to variations in seeing conditions. 



060607 A, 07 1 03 1 , 0803 1 0, 0803 1 9B and 0804 1 3 A) correspond- 
ing to Az = 9.4 that are not in the Prochter's sample. Therefore 
more than two third of the redshift path is new. In addition, the 
UVES sample is homogeneous (similar resolution, same instru- 
ment, similar signal-to-noise ratio) and more than doubles the 
sample used by Sudilovsky et al. (20071, which is included in 
our sample, for the Mg n statistics. 

The second sample we consider (the overall sample) is 
formed by adding observations from the literature (see Table [4]i 
to the UVES sample (see Section 3.2). The sample gathers ob- 
servations of 26 GRBs for Az = 31.55. It therefore doubles the 
statistics of |Prochter et al.| ( |2006b"l > 



In the f ollowing we use solar abundances from Gre vesse 

|etaT1 ( |2007l l. 

3. Number density of Mg n absorbers 

3.1. The UVES sample 

For each line of sight we searched by eye the spectrum for Mg n 
absorbers outside the Lyman-Q- forest considering all Mg n com- 
ponents within 500 km s _1 as a single system. Table pjsumma- 
rizes the results. Columns 1 to 8 give, respectively, the name of 
the GRB, its redshift, the redshiftpaths along the line of sight 
for W r ji m > 0.3 and 1 A (see Eq.njl, the redshift of the Mg n ab- 
sorber, the rest equivalent width ofthe Mg H/12796 transition and 
the velocity relative to the GRB redshift. The last column gives 
comments on peculiar systems if need be. When a Mg n/12796 
line is blended either with a sky feature or an absorption from an- 
other system, we fit simultaneously Voigt profile components to 
the Mg n doublet and the contaminating absorption and derive 
characteristics of the Mg n/12796 absorption from the fit. The 
sky spectrum at the position of a feature is obtained from other 
UVES spectra in which the metal absorption is not present. The 
upper limits include also the contaminating absorption. 

The rest equivalent width detection limit (at any given sta- 
tistical level) is calculated at each redshift along the spectrum 
using the following equation (Tytle r & Fan|19 94): 



W, 



r,lim 



UM^ 5 
SNR 



AA (A), 



(1) 



where Ml is the number of pixels the line is expected to cover, 
U is the number of rms-intervals (or <x) defining the statistical 
significance of the detection limit, SNR is the signal-to-noise 
ratio at the corresponding wavelength and AA is the FWHM of 
the resolution element. We will apply this detection limit to the 
Mg n/12796 transition. Using Ml = 5 and U = 6 (so that each 
transition of the Mg n doublet is detected at more than 3 cr), it can 
be seen that for a typical signal-to-noise ratio of SNR ~ 10 and 
a typical FWHM of AA = 0.13 A (R = 43000 and A = 5600 A), 
the formula gives a detection limit of ~0.17 A which is much 
smaller than the equivalent width limit we will use in our statis- 
tical studies below (1 and 0.3 A). 

We then can compute the redshift path over which a line of 
a given equivalent width would be detected in our data. |Prochter| 
et al. ( 2006b i hmited their analysis to the redshift range 0.359-2 
and compared their results to those of the SDSS QSO survey 
reported by |Prochter et al.| ( |2006a| ) However, the SDS S QSO 
Mgll survey extends to z max ~ 2.3 ( Nestor et al.||2005 l so that 
there is no reason to limit our analysis to z = 2 and we will use 
the same redshift limits as Nestor et al. (2005), Zstart = 

0.366 

and z m ax = 2.27 instead. As in all QSO surveys and follow- 
ing Prochter et al. (2006a), we will exclude along each line- 
of-sight the redshift range within an ejection velocity of 3000 
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GRB 


ZGRB 


Az 

W r >0.3A 


Az 
W r >1A 


■Sabs 


W,. 0*2796) 
(A) 


v CJ 
(km/s) 


Remarks 




021004 


2.3295 


1.754 


1.756 


0.5550 
1.3800 
1.6026 


0.248 ± 0.025 
1.637 ±0.020 
1.407 ± 0.024 


~ 192000 
~ 97000 
~ 72000 


blended with A1II1670 at z = 


1.6026 


050730 


3.9687 


1.278 


1.298 


1.7732 
2.2531 


0.927 ± 0.030 
< 0.783 (0.650) 


~ 157000 
~ 120000 


sky contamination subtracted 




050820A 


2.6147 


1.843 


1.845 


0.6896 
0.6915 
1.4288 
1.6204 


0.089 ± 0.007 
2.874 ± 0.007 
1.323 ±0.023 
0.277 ± 0.024 


~ 192000 
~ 192000 
~ 113000 
~ 93000 














2.3598 


< 0.424(0.306) 


~ 22000 


contribution by FeII2600 at z 
subtracted 


= 2.6147 


050922C 


2.1996 


1.679 


1.682 


0.6369 
1.1076 
1.5670 


0.179 ± 0.018 
0.677 ± 0.029 
< 0.102 (0.08) 


~ 175000 
~ 118000 
~ 62000 


sky contamination subtracted 




060418 


1.4900 


1.242 


1.265 


0.6026 
0.6559 
1.1070 


1.293 ± 0.010 
1.033 ±0.006 
1.844 ± 0.014 


~ 124000 
~ 116000 
~ 50000 






060607A 


3.0748 


1.710 


1.713 


1.5103 
1.8033 
2.2783 


0.204 ±0.011 
1.854 ±0.006 
0.343 ± 0.058 


~ 135000 
~ 107000 
~ 64000 






071031 


2.6922 


1.727 


1.789 


1.0743 
1.9520 


0.330 ±0.016 

\j.\yyL + U.U1 + 

0.804 ±0.016 


~ 156000 
~ Q7nnn 

~ 66000 






080310 


2.4272 


1.841 


1.841 


1.1788 
1.6711 


0.047 ± 0.024 
0.421 ± 0.012 


~ 127000 
~ 73000 






0803 19B 


0.9378 


0.790 


0.790 


0.5308 
0.5662 
0.7154 
0.7608 


0.614 ±0.001 
0.083 ± 0.003 
1.482 ±0.001 
0.108 ±0.002 


~ 69000 
~ 63000 
~ 36000 
~ 29000 






080413A 


2.4346 


1.583 


1.650 


2.1210 


0.768 ± 0.259 


~ 29000 







km s" 1 from the GRB redshift. Table [5] lists the mean redshift, 
(Zabs), the total redshift paths obtained considering these red- 
shift limits, and the number of systems detected in our sam- 
ple over these redshift paths for different W T limits or ranges: 
W t > 0.3 A ( 'all systems '), W r > 1 .0 A ( 'strong systems ') and 
0.3 < W r < 1.0 A('weak systems'). The total redshift paths 
for the all and strong samples are Az = 13.79 and 13.94, re- 
spectively, for an observed number, N , , of 18 and 9 sys- 
tems detected, corresponding to redshift number densities of 
dn/dz = 1.31 ± 0.31, 0.65+0.22 for all and strong Mg n sys- 
tems, respectively. 



Table 3. Number of Mg n systems and redshift paths 



W r (/!2796) 


>0.3A 


>1 A 


> 0.3 and <1 A 


(Zabs) 


1.34 


1.11 


1.57 


Redshift path 


13.79 


13.94 


13.79 


N 1f [ (UVES sample) 


18 


9 


9 


jyMgll 
iv exp 


11.98 (±3.46) 


4.83 (±2.20) 


7.21 (±2.68) 


(Nestor et al., 2005) 








iv exp 




4.00 (±2.00) 




(Prochter et al., 2006) 









Fig. 1. Redshift path density g(z) of the UVES sample for W r > 
0.3 A (black dashed line) and W t > 1 .0 A (red line). 



In Table[3] we also report the number of Mg n absorbers that 
would be expected along lines of sight toward QSOs over the 
same redshift path, A^jf 1 . To calculate these numbers, we deter- 
mine the redshift path density, g(z), of the UVES GRB sample 
(see Fig. [T| and combine it with the density of QSO absorbers 
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Mgll systems, W r >0.3A 



QSOs (Nestor et 


al. 2005) 


GRBs.UVHSsa 























Mgll systems, W r > 1.0A 




Fig.|2]and reported in Table[3] Errors are assumed to be poisso- 



nian and are scaled as 



Mgll 



exp 



As a verification, we note from Fig. 13 of Nestor et al. ( 2005| l 
that the number density, dn/dz, of W T > 0.3 A systems is almost 
independent of the redshift. Multiplying the dn/dz = 0.783 value 
reported there by the UVES 'all systems' Az = 13.79 we obtain 
an expected total number of 9.78, in agreement with the value 
reported in Tab. [3] for systems along QSO los (i.e. 10.54). The 
same test performed using the number density values found by 
Steide l & Sargent ( 1992) give consistent results. 



Fig.|2|and Table 3 show that the excess of strong Mg n sys- 
tems along GRB los compared to QSO los is significant at more 
than 2cr (slightly less than 2<x for the strong Mg n systems if the 
[Nestor et al. 2005 function is used), but it is more than a fac- 
tor of ~ 2 lower than what is found by Prochter et al. (2006b). 
The number of weak systems is consistent within ~ 1 <x to that 
expected for QSO los. 

3.2. The overall sample 



Mgll systems, 0.3 A < W r < 1 .0 A 



Mgll systems, W r > 1.0 A, 'overall sample' 



QSOs (Nestor et 


al,2005) 




CRBs.UVESsa 





























Fig. 2. Comparison between the cumulative distribution of Mg n 
systems detected along the UVES GRB los (red) and the one 
expected along QSO los following Eq. 2 (black solid) or Eq. 3 
(black dashed), for: all (top panel), strong (middle panel) and 
weak (bottom panel) systems. 



. QSOs (Nestor et al. 2005) 
QSOs (Prochter et al. 2006) 
GRBs, overall sample 




Fig. 3. Comparison between the cumulative distribution of 
strong Mg n systems in the overall sample of GRB los (red) and 
the one expected along QSO los following Eq. 2 (black solid) or 
Eq. 3 (black dashed). 



per unit redshift as observed by Nestor et al. ( 2005 \ and Prochter 
et al. ( 2006a i (for the strong systems only) in the SDSS survey: 



n 



dz. 



(2) 



Both |Nestor et al.|P005) and |Prochter et al.| ( |2006al > showed 
that the traditional parametrization of dn/dz as a simple pow- 
erlaw n (l + z) y does not provide a good fit to the SDSS data. 
Therefore we will use their empirical fits to the redshift number 
density. [Nestor et aL] ( |2005[ > give 



dn 

dz 



I 



dn 
dWt 



(z)dW : 



f n\z) 
J W*(z)' 



w,/w 



dW, 



(3) 



where both n* and W* vary with redshift as power laws: n* = 
1.001+0.132(l+z) Q - 22 6±Q170 , W* = 0.443+0.032(1 +z)°- 634±0097 . 
|Prochter et al.] ( 2006b) derive for the strong Mgll systems: 



dn/dz = -0.026 + 0.374z - 0.145z 2 + 0.026z 3 



(4) 



The results of the calculations for the different W r ranges 
(Wr > 0.3 A, Wr > 1.0 A and 0.3 < W r < 1.0 A) are shown in 



In order to increase the statistics of strong Mg n systems, we 
added to the UVES sample both high and low resolution GRB af- 
terglow spectra published in the literature. The resulting sample 
is composed of 26 los (see Table[4]for details on los that are not 
part of the UVES sample). Since this sample includes many low 
resolution spectra, we study the statistics of strong (W r > 1 A) 
Mg ii absorbers only, using the same redshift limits as in the pre- 
vious Section. Three lines of sight of this sample (GRB 991216, 
GRB 00 0926 and GRB 030429) were not used by |Prochter etaL 
(2006b) because the low spectral resolution of these spectra does 
not allow to resolve the Mg n doublet. However, such a low reso- 
lution does not prevent the detection of strong systems, although 
the doublet is blended. In addition, the total equivalent width 
of the doublets detected along these los is larger than 3 A so 
that we are confident that W t jn% is larger than 1 A. In any case, 
as detection of Mg n systems is more difficult at lower resolu- 
tion, including these lines of sight could only underestimate their 
number density. 

The total number of strong Mg n systems is N = 22 and 
the redshift path is Az = 31.55. This leads to a number density 
dn/dz = 0.70 + 0.15. We use the g(z) function of this enlarged 
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Table 4. GRB los available from the literature. 



5 



GRB 


ZGRB 


^start 


<-tMld 


Zabs 


tV r (2796) 
(A) 


Av ej 

(km/s) 


Reference 


991216 


1.022 


0.366 


1.002 


0.770° 


4.0 ± 0.8" 


~ 40000 


3 










0.803 


6.1 ±0.7* 


~ 34000 




000926 


2.038 


0.616 


2.008 








1 


010222 


1.477 


0.430 


1.452 


0.927 


1.00 ±0.14 


~ 74000 


2 










1.156 


2.49 ± 0.08 


~ 41000 




011211 


2.142 


0.366 


1.932 








3 


020405 


0.695 


0.366 


0.678 


0.472 


1.1 ±0.3 


~ 42000 


4 


020813 


1.255 


0.366 


1.232 


1.224 


1.67 ±0.02 


-4000 


5 


030226 


1.986 


0.366 


1.956 








6 


030323 


3.372 


0.824 


1.646 








7 


030328 


1.522 


0.366 


1.497 








14 


030429 


2.66 


0.620 


1.861 


0.8418 


3.3 ±0.4* 


~ 179000 


15 


050505 


4.275 


1.414 


2.27 


1.695 


1.98 


~ 176000 


8 










2.265 


1.74 


~ 134000 




050908 


3.35 


0.814 


2.27 


1.548 


1.336 ±0.107 


~ 147000 


9 


051111 


1.55 


0.488 


1.524 


1.190 


1.599 ±0.007 


~ 45000 


9,13 


060206 


4.048 


1.210 


2.27 


2.26 


1.60 


~ 123000 


10 


060526 


3.221 


0.836 


2.27 








11 


071003 


1.604 


0.366 


1.578 


0.372 


2.48 ± 0.20 


~ 170000 


12 



1 



Castro et al.]{2003f ; 2:|Mirabal et aL 



(20021; 3 



6 :|Klose et al.l(|2604b ; 7: |Vreeswijk et al 



2004 1; 8: 



ll: |ThoneetaTTp008) ; 12: Perley et al.|( |2008| >; 13:[T 



Vreeswijk et al-H2006| >; 
Bergeretal.|(|2006^;9:| 



5:|Barth et al.|d2003|»; 4:|Masetti et al.lt 2003 1 



ao et al.| ( |2067[ >; 14: Maiorano et al. (2006); 15: Jakobsson et al. (2004). 
"The very low resolution of the GRB 991216 spectrum makes the z=0.77 absorber identification uncertain. 6 The W T values refer to the total 

equivalent width of the Mgll doublet. 



sample to compute the total number of strong systems expected 
for a simi lar QSO sample following the same method as used in 
Sect. |3~ij We find N™f = 10.41 + 3.23 andA^ 11 = 8.62 ±2.94, 
using Eq. 3 and 4, respectively, that is 2.1 ± 0.6 and 2.6 + 0.8 
times less than the number found along GRB los (see Fig. [3]). 

The excess of strong Mg n systems along GRB lines of sight 
for this enlarged sample is confirmed at a ~ 3<x statistical sig- 
nificance. The excess found is higher than for the UVES sample 
but still a factor of ~ 2 lower than what was previously reported 
by Prochter et al. ( 2006b i from a sample with a smaller red- 



shift path. The number densities resulting from the two studies 
are different by no more than 2cr. For consistency we also per- 
formed our analysis considering only the smaller sample used 
by Prochter et al. ( 2006b i. In this case, the results obtained are 
similar to those found by these authors. The redshift path of our 
overall sample is twice as large as that used by [Prochter et al. 
(2006b), therefore the factor of 2 excess we find in this study 
has a higher statistical significance. 



3.3. Number density evolution 

We divided both the UVES and the overall sample in three red- 
shift bins and calculated dn/dz for each bin. Fig. |4] shows the 
number of systems per redshift bin. While the total number of 
systems and the number of weak systems have a comparable 
redshift evolution in GRB and QSO lines of sight, the strong 
systems happen to have a different evolution in GRB los. The 
excess of strong systems in GRBs is particularly pronounced at 
low redshift, up to z ~ 1.6. We performed a KS test for each of 
the three cases to assess the similarity of the redshift distribu- 
tion of Mg ii systems along GRB and QSO los. There is a 90.7% 
and 23.5% chance that the weak and all Mg n absorber sam- 
ples in QSOs and GRBs are drawn from the same population. 
The probability for the strong systems is 20.1% for the UVES 



sample, but it decreases to ~ 2% when considering the overall 
sample. 

This apparent excess of strong systems in the low-redshift 
bin could indicate that some amplification bias due to lensing is 
at work. Indeed the effect of lensing should be larger in case 
the deflecting mass is at smaller redshift. The lensing optical 
depth for GRBs at redshift zgrb >> 1 and zgrb < 1 is maxi- 
mal for a lens at respectively z\ ~ 0.7 and zgrb/2. However, we 
find that only 47% of GRBs with strong foreground absorbers 
in our sample have at least one strong Mg n system located 
in the range within which the optical depth decreases to about 
half its maximum value ( Sudilovsky et al.||2007| [Turner| [l980 



Sme tte et al.|1997 l. This indicates that if amplification by lens 
ing is the correct explanation, the effect must be weak and subtle. 
Porciani et al. ( 2007) used the optical afterglow luminosities re- 
ported by Nardini et al. (2006 ) to show that GRB afterglows with 
more than one absorbers are brigther than others by a factor 1.7. 
However this correlation has not been detected by Sudilovsky 
et al. ( 2QQ7\ when using the afterglow B-band absolute magni- 



tudes obtained by |Kann etal| (2006 2007]). 

More data and larger samples are obviously needed to con- 
clude on this issue. 

4. The population of strong Mg n systems 

Both the results on the UVES homogeneous sample and those 
obtained using the overall sample confirm, although to a smaller 
extent, the excess of strong Mg n absorbers along GRB los first 
reported by |Prochter et al.| ([2006b). To understand the reason of 
the discrepancy in the number density of "strong" systems it is 
therefore important to study in more details these systems and 
to derive their physical characteristics. This is possible using the 
UVES data. The main question we would like to address here is 
whether there is any reason to believe that GRB and QSO strong 
absorbers are not drawn from the same population. 
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Fig. 4. Number density evolution of Mg n systems detected 
along GRB los for: first panel: all, second panel: strong, third 
panel: weak systems of the UVES sample; bottom panel: strong 
systems of the overall sample. The solid line represents the evo- 
lution of the Mgll number density derived from the Mg n survey 
in the SDSS QSO by |Nestor et aLlpOOB) (see Eq. 3). 



4. 1 . Equivalent width distribution 

We show in Fig. [5] the comparison between the normalized W r 
distribution of all Mg n systems with W T > 0.3 A detected in the 
UVES sample and the one reported by Nestor et al. ( 2005 1 for 
the Mg n systems along the QSO los in the SDSS survey. 

The KS tests give a 27.1% chance that the GRB and QSO 
distributions are drawn from the same population for the UVES 
sample. This result extends to lower W r values the conclusion 
by Porciani et al.| (|2007 ) that the two distributions are similar, 
arguing against the idea that the excess of Mg n systems could 
be related to the internal structure of the intervening clouds. 




Fig. 5. W r distribution of the Mg n systems with W T > 0.3 A de- 
tected in the UVES sample. The solid curve represents the cor- 



responding distribution for QSO Mgll absorbers (Nestor et al. 
2005 ). Both distributions have been normalized such that their 
underlying areas are equal to 1 . 
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Fig. 6. The W r > 1 .0 A Mg n absorption systems detected in the 
UVES spectra, the Mg n ,12796 (top) and /12803 (bottom) ab- 
sorption lines are displayed for each system on a relative veloc- 
ity scale with km s _I centered at the redshift reported in Table 
1. 



4.2. The velocity spread of strong Mg n systems 

The profiles of the two transitions of the W l > 1 A Mg n systems 
found in the UVES spectra are plotted on a velocity scale in 
Fig. [6] The profiles are complex, spread over at least 200 km s _1 
and up to ~ 600 km s _I and show a highly clumpy structure. 

We measure the velocity spread of the Mg n systems with 
W r > 1.0 A in the UVES sample following [Ledoux et al (2006 1. 



We therefore use a moderately saturated low ionization absorp- 
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Fig. 7. The solid line histogram represents the distribution of the 
velocities of the strong Mg n systems in the UVES sample cal- 
culated following |Ledoux et al. ( 2006| l. The thick dashed his- 



togram refers to the 27 SDSS Mg n systems observed at high 
spectral resolution ( [Prochter e t al. 2006a]!. 



tion line (e.g. Fell/12344, Sill/11526) and measure the velocity 
difference between the points of the absorption profiles at which 
5% and 95% of the absorption occurs. This method is defined 
in order to measure the velocity width of the bulk of gas, avoid- 
ing contamination by satellite components which have negligi- 
ble contribution to the total metal column density. Note that us- 
ing this definition implies that the measured velocity spread is 
usually smaller than the spread of the Mg n profile which is 
strongly saturated. In case no moderately saturated absorption 
line is available, we use the mean value of the velocity widths 
calculated both from a saturated line and an optically thin line. 
Results are given in column 9 of Table [6] 

The velocity-spread distribution is shown for UVES systems 
with W r > 1 A in Fig. [^together with the SDSS QSO distribution 
from 27 SDSS systems with W r > 1 A observed at high spec- 
tral resolution ( [Prochter et al.|2006a[ >. Although the statistics are 
small, the UVES and SDSS distributions are statistically simi- 
lar. A KS test gives 65% chance that the two samples are drawn 
from the same population (see also Cucchi ara et al.|2 008 ). 



4.3. Dust extinction 

It has been proposed that a dust bias could possibly affect the 
statistics of strong Mg n systems. Indeed, if part of the popu- 
lation of strong Mg n systems contains a substantial amount of 
dust, then the corresponding lines of sight could be missed in 
QSO surveys because of the attenuation of the quasar whereas 
GRBs being intrinsically brighter, the same lines of sight are 
not missed when observing GRBs. The existence of a dust ex- 
tinction bias in DLA surveys is still a debated topic among the 
QSO community although observations of radio selected QSO 
los ( |Ellison eT al. 2004) seem to show that, if any, this effect is 
probably small (see also Pontzen & Pettini 2008 ). 

We can use our UVES lines of sight to estimate the dust con- 
tent of strong Mg n systems from the depletion of iron com- 
pared to other non-depleted species as is usually done in DLAs. 
We have therefore searched for both Fe n and Zn n absorption 
lines (or Si n, in case the Zn n lines are not available) associ- 
ated to the strongest Mg n systems present in the UVES spec- 
tra. Because the spectra do not always cover the relevant wave- 



length range, the associated column densities could be measured 
only for 4 out of the 10 systems (see Table [5j. We estimate the 
depletion factor, and therefore the presence of dust, from the 
metallicity ratio of iron (a species heavily depleted into dust 
grains in the ISM of our Galaxy) to zinc (that is little depleted), 
[Fe/Zn] = (Fe/Zn)-(Fe/Zn)0 (or [Fe/Si]). We also determine the 
iron dust phase column density (N^ st ) using the formula given 



by |Vladilo et al.| ([2006 ) and from this we infer an upper limit 
on the corresponding flux attenuation Ay from their Fig. 4. We 



used also the estimate given by (Bohlin et al. 1978 Prochaska & 
|Wolfe|2002| i: 



10* 



10 21 x 10 [X/H] Q 



(1 - 10 [Fe/xl ) 



(5) 



with k = 10 [X/H] (1 - 10 [Fe/x] ) representing the dust-to-gas 
ratio and X corresponding to Zn or Si if Zn is not available. 

Results are reported in Table [5] It can be seen that al- 
though depletion of iron can be significant, the corresponding 
attenuation is modest because the column densities of metal- 
lic species are relatively small owing to low metallicities (see 
also Table 5). Indeed we find high dust depletion in two sys- 
tems at z a bs = 1.107 toward GRB 060418 (see also |Ellison et al.| 
[20061 ) and at z abs = 1.6026 toward GRB 021004. Their col- 
umn densities are low however so that the inferred A v does 



not exceed values typically found for QSO los (see Vladilo 
et~aI1l2006l and |Prantzos & Boissier||20T)0"] ). |Cucchiara eTal 
(2008|) donot detect the Znll absorption lines for the systems 



at Zabs = 1.6026 toward GRB 021004. They report an upper limit 
of W r < 0.016 A whereas we find W r = 0.07 A. The detection 
of Zn II(/1/12026,2062) for the central component seems robust, 
as it is shown in Fig. [8] In addition, the Cr II(/L12056,2062) lines 
are also possibly detected. 

We note that for GRB 021004 and GRB 060607 A we could 
estimate the attenuation for ALL strong Mg n systems identified 
along the los. Therefore our results do not show evidences to 
support the idea that a bias due to the presence of dust in strong 
Mg ii systems could be an explanation for the overabundance of 
the strong Mg n absorbers along GRB los. On the other hand, the 
fact that 2 out of 4 of the selected systems have a dust depletion 
higher than the values usually found for QSO DLAs (Meiring 
|et al.|p006) ) supports the fact that these systems probably arise 
in the central part of massive halos where the probability to find 
cold gas is expected to be higher. 



5. Estimating the H i column density 

A key parameter to characterize an absorber is the corresponding 
H i column density. Unfortunately, the H i Lyman-a absorption 
line of most of the systems is located below the atmospheric 
cut-off and is unobservable from the ground. If we want to char- 
acterize the systems with their H i column density or at least an 
estimate of it, we have to infer it indirectly. For this we will as- 
sume that the strong systems seen in front of GRBs and QSOs 
are cosmological and drawn from the same population. We be- 
lieve that the results discussed in the previous section justify this 
assumption. 

We use the velocity-metallicity correlation found by |Ledoux| 



et al. (2006 1 to estimate the expected metallicity of the systems 



in the UVES sample with W T > 1.0 A, assuming that the corre- 
lation is valid also for sub-DLA systems (Peroux et al. 2008). 
The relation was linearly extrapolated to the average redshift of 
the sample of systems with W t > 1.0 A , (z) = 1.11, giving 
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Table 5. Iron to Zinc or Silicon ratio and extinction estimate for 4 strong Mg n systems. 



A'(Fell) 

(cm" 2 ) 



A'(Znll) 

(cm- 2 ) 



Af(Sill) 

(cm- 2 ) 



[Fe/Zn] 



[Fe/Si] 



N£ s< A v 
(cm" 2 ) 



GRB021004 1.3800 15.09 ± 0.05 



15.19 ±0.09 



-0.05 ±0.10 14.18 ±0.23 



GRB021004 1.6026 14.60 ± 0.03 12.78 ± 0.02 



-1.03 ±0.04 



15.59 ± 0.02 < 0.2 



GRB060418 1.1070 14.69 ±0.01 12.87 ± 0.03 



-1.05 ±0.03 



15.68 ± 0.03 < 0.2 



GRB060607A 1.8033 14.07 ± 0.03 



14.36 ±0.10 



-0.24 ±0.10 13.94 ±0.17 ~0 



Table 6. Characteristics of Mg n systems with W r > 1 .0 A in the UVES sample. 



GRB 


<-abs 




yyMgIA2S52 




A^znn 


AW 


AW 


"Av 


b^ v MgIlA2796 






(A) 


(A) 


(A) 


-2 

cm 


cm -2 


-0 

cm - 


km s" 1 


km s 1 


021004 


1.380 


1.637 


0.190 


0.972 




15.19 ±0.09 


15.09 ±0.05 


170 


89 


021004 


1.603 


1.407 


0.366 


0.737 


12.78 ± 0.02 




14.60 ±0.03 


164 


51 


050820A 


0.692 


2.874 


N/A c 


N/A c 






N/A c 




228 


050820A 


1.429 


1.323 


0.488 


0.601 






14.34 ± 0.03 


271 


98 


060418 


0.603 


1.293 


0.361 


0.989 






16.43 ±0.04 rf 


76 


130 


060418 


0.656 


1.033 


0.078 


0.486 






13.96 ±0.02 


136 


123 


060418 


1.107 


1.844 


0.483 


1.080 


12.87 ±0.03 




14.69 ± 0.01 


119 


221 


060607A 


1.803 


1.854 


> 0.226 


0.825 




14.36 ±0.1 


14.07 ± 0.03 


333 


90 


0803 19B 


0.715 


1.482 


0.303 


0.697 






14.00 ± 0.01 


354 


142 



" Velocities are measured following Ledoux et al.| ( |2006) ; b Velocities are measured using the MgII/12796 absorption considering only the central 
part of the system, following the method recommended by |Ellison et al.| ( p009| l; c Lines redshifted on top of the Ly-a absorption associated to the 
GRB or in the Ly-a forest; d Saturated line. 



[X/H]= 1.431ogAv - 3.78. We infer the hydrogen column densi- 
ties dividing the zinc, silicon or iron column densities measured 
in the UVES spectrum by the metallicity. Note that in most of 
the cases only Fe n available (see Table 5). The A^(H i) column 
density derived in these cases should be considered a lower limit 
because iron can be depleted onto dust-grains. The results are 
shown in Table [7] (columns 3 and 4). The error on the proce- 
dure should be of the order of 0.5 dex (see |Ledoux et al.|2 006 ). 
We insist on the fact that our aim is not to derive an exact H i 
column density for each system but rather to estimate the over- 
all nature of the systems. We find that at least 3 of the 9 sys- 
tems with W r > 1 A could be DLAs. In any case, and even 
if we consider that we systematically overestimate the column 
density by 0.5 dex, a large fraction of the systems should have 
logAMTIi) > 19. 

Another method to infer the presence of DLAs among Mg n 
absorbers is to use the 'D-index' (Ell ison et al.||2009] l. We cal- 
culate the D-index for the systems in our sample following the 
recommended D-index definition by Ellison et al. ( 2009) 1 and us- 
ing the formula: 



(column 5). Note that if we do not include the iron column den- 
sity term in the D-index calculation, we obtain similar results. 

We find that all the 9 systems have D > 7. [Ellison et al. 



(2009) find that 57+J^% of the systems having D > 7 are DLAs. 



This results applied to our sample imply that at least 5 systems 
are DLAs. We can compare the number of DLAs per unit red- 
shift («dla) found for the QSO los by |Rao et aL] ( |2006| l to that of 
our GRB sample. «dla is defined as the product of the Mg n sys- 
tems number density and the fraction of DLAs in a Mg n sample. 
IRao et al.|(|2006|l found «dla,qso ~ 0.1 (with errors of the order 
of 0.02). If we assume that three to five systems in our sample 
are DLAs we find «dla,grb ~ 0.22 - 0.36 . This means that the 
number of DLAs is at least two times larger along GRB lines of 
sight as compared to QSO lines of sight. 

The above estimate of log A^(H i) can be considered as highly 
uncertain and the identification of a few of the systems as DLAs 
can be questioned. 



D = 



W r (MgIU2796)(A) 
Av(MgIL12796)(km/s) 



X x 1000. 



(6) 



The width of the central part of the Mg n/12796 absorption 
system, Av(MgILl2796) in km/s, is reported in column 10 of 
Table|6]while the resulting D-index values are reported in Table|7] 



All this seems to indicate that GRBs favor lines of sight with 
an excess of DLAs. Since these systems are more likely to be 
located in the central parts of massive halos, this may again favor 
the idea that there exists a bias in GRB observations towards 
GRBs with afterglows brighter because they are subject to some 
lensing amplification. 
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Fig. 8. Fe 11(^1608,2344,2374), Zn 11(^2026,2062) and 
Crll(/l/t2056,2062) transition lines at z a t, s = 1.6026 in the spec- 
tra of the afterglow of GRB 021004, shown in velocity space. 
The dashed line correspond to z a b s = 1.6026. The Znll absorp- 
tion lines are clearly detected for the central component of the 
system. 



Table 7. Inferred metallicity, Nm and D-index of Mg n systems 
with W r > 1.0 A in the UVES sample. 



GRB 


Zabs 


[X/Hf 




D-index' 


021004 


1.380 


-0.63 


20.3 


8.32 + 0.22 


021004 


1.603 


-0.65 


20.8 


9.99 ± 0.47 


050820A 


0.692 






9.0 ± 0.03 rf 


050820A 


1.429 


-0.34 


> 19.2 


8.13 ±0.22 


060418 


0.603 


-1.13 


> 22.0 


10.39 ± 0.09 f 


060418 


0.656 


-0.77 


> 19.2 


7.51 ±0.05 


060418 


1.107 


-0.85 


21.1 


8.17 ±0.06 


060607A 


1.803 


-0.21 


19.0 


9.12 ±0.06 


0803 19B 


0.715 


-0.17 


> 18.7 


8.06 ±0.01 



Metallicities infer red using the velocity-metallicity relation found by 
Ledoux et al.|{2006|; b Nm v alues inferred using the velocity-metallicity 
relation found by |Ledoux et al.| ([2006 ). When only the iron column den- 
sity is available, the N m value has to be considered has a lower limit due 
to possible dust extinction; c D-index calculated using Eq.[6] following 
the method recommended by Ellison et al. 1 2009 ; d Z)-index calculated 
without including the iron column density term (see |Ellison et al.|2 009 1; 
e iron lines are saturated. The D-'mdex calculated without including the 
iron column density term would be D = 9.48 ± 0.08; 



6. Observed sub-DLA absorbers 

The UVES spectra often cover a substantial part of the Lyman- 
a forest in front of the GRB. It is therefore possible to search 
directly for strong H i Lyman-a absorption lines correspond- 
ing to DLAs or, more generally, sub-DLAs. This can only 
be performed for 8 of the 10 GRBs in the sample: the red- 
shifts of GRB 060418 and GRB 0803 19B are unfortunately too 
low to allow the detection of the Lyman-ff absorptions in the 



UVES spectral range. We exclude from our search the DLA at 
the GRB redshift, which is believed to be associated with the 
close surrounding of the GRB. We find additional (sub-)DLAs 
(Fig.[9]) along the los of GRB 050730 (s ee also|Chen et al.|2005), 
GRB 050820, GRB 050922C (see also |Piranomonte et al.|2 008 ) 
and GRB 060607 A (Fig. 10 1. The former three systems are sim- 



ple, with a single Lyman-a- component. The system detected at 
2.9374 toward GRB 060607A is more complex and cannot be fit 
with a single component (sub-) DLA. The profile is made of two 
main clumps at z — 2.9322 and 2.9374. The latter is a blend of 
the four O i components (see top panel in Fig. 9). The constraints 
on the H i column densities come mostly from the red wing of 
the Lyman-Q' line and the structures seen in the Lyman-y line. 
We find that the two main H i components at z = 2.9374 and 
2.9322 have log N(H i) = 19.4 and 19.0 ± 0.1, respectively. It is 
worth noting that the limit on the O i component at z = 2.9322 is 
log A^(0 1) < 12.6 implying that the metallicity in the cloud could 
be as low as [O/H] < -2.9 which would be the lowest metallic- 
ity ever observed yet in such a system. Some ionization correc- 
tion could be necessary however in case the oxygen equilibrium 
is displaced toward high excitation species by a hard ionizing 
spectrum. 

Table[8] shows the column densities and abundances of the 
systems assuming that all elements are in their neutral or singly 
ionized state. The values are similar to those usually found for 
(sub)-DLAs along QSO los (see for example Peroux et al. 2003 ). 

The presence of these systems allows us to extend the search 
for strong systems along GRB lines of sight to higher redshift. 
The redshift paths (see Table|8]l have been calculated considering 
Zstart as the redshift of an absorber, for which the corresponding 
Lyman-Q' line would be redshifted to the same wavelength as the 
Lyman-/? line of the GRB afterglow. We therefore avoid confu- 
sion with the Lyman-/? forest. z en d is fixed at 3000 km/s from the 
GRB redshift. The only exception to this rule is GRB 050730 for 
which there is a gap in the spectrum located at about 3000 km/s 
from the GRB redshift and starting at A = 5096 A. We use this 
wavelength to fix z max . The corresponding total redshift path is 
Az = 4.17. The sub-DLAs systems at z = 2.9322 and 2.9374 
along GRB 060607 A are separated by ~400 km s , therefore 
for the statistical study we consider them as a single system. The 
resulting number density for systems with log A^(H i)> 19.0 is 
0.96 ± 0.48 for an average redshift of (z) = 2.08. At this redshift 
a value of about 0.5 + 0.2 is expected from QSO los (|Peroux 



et al.|2005) . There could be therefore an excess of such systems 
but the statistics is obviously poor. 

It is intriguing however that the four systems we detect are all 
in the half of the Lyman-a- forest closest to the GRB, with ejec- 
tion velocities of about 25000, 22000, 12000 and 1 1000 km s" 1 
for GRB 050730, 050820A, 050922C and 060607A respec- 
tively. The corresponding ejection velocities found for the strong 
Mg ii systems (see column 7 of Table 2) are larger than 
35000 km/s, making an origin local to the GRB unlikely for the 
strong Mg n systems. This fact therefore does not favor the pos- 
sibility that the excess of strong Mg n absorbers along GRB los 
is due to some ejected material present in the GRB environment. 
The much lower ejection velocities found for the sub-DLA ab- 
sorbers may indicate that for these systems the excess is not of 
the same origin as the excess of Mg n systems and that part of 
the gas in these systems could have been ejected by the GRB. 

7. Conclusion 

We have taken advantage of UVES observations of GRB after- 
glows obtained over the past six years to build an homogeneous 
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Fig. 9. Foreground (sub)-DLAs. Top panel: Sub-damped Lyman- 
a absorption at z = 3.5655 toward GRB 050730. Middle 
panel: Damped Lyman-or absorption at z = 2.3598 toward 
GRB 050820. Bottom panel: Sub-damped Lyman-ff absorption 
at z = 2.0778 toward GRB 050922C (this system has been fitted 
together with the DLA associated to the GRB at z = 2.1999). 
The red area correspond to the fit result covering the column 
density error range. 



sample of lines of sight observed at high spectral resolution and 
good signal-to-noise ratio (SNR>10). We used these observa- 
tions to re-investigate the claimed excess of Mg n systems along 
GRB lines of sight extending the study to smaller equivalent 
widths. We also used these data to derive intrinsic physical prop- 
erties of these systems. 

Considering the redshift ranges 0.37 < z < 2.27 of the SDSS 
survey used for the QSO statistics, we find an excess of strong 
intervening Mgll systems along the 10 GRB lines of sight ob- 
served by UVES of a factor of ~ 2 compared to QSO lines of 
sigth. This excess is significant at ~ 2cr. Thanks to the quality of 
the UVES data it has been possible also to consider the statistics 
of the weak absorbers (0.3 < W T < 1.0 A). We find that the num- 
ber of weak absorbers is similar along GRBs and QSOs lines of 
sight. 

We increase the absorption length for strong systems to 
Az = 31.5 using intermediate and low resolution observations 
reported in the literature. The excess of strong Mgll interven- 
ing systems of a factor of ~ 2 is confirmed at 3<x significance. 
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Fig. 10. The sub-DLA system at z = 2.9374 along the los to- 
ward GRB 060607 A. The top panel shows the fit of the O l11302 
absorption profile, while the following panels report the fit re- 
sults of the Lyman-y, Lyman-/5 and Lyman-*? absorption lines 
(from top to bottom, respectively). The red lines show the main 
z = 2.9374 component (different lines include errors) mostly 
constrained by the red wing of the Lyman-o' absorption. In 
blue is shown the total fit including the sub-DLA component 
at z = 2.9322. 



We therefore strenghten the evidence that the number of strong 
sytems is larger along GRB lines of sight, even if this excess is 
less than what has been claimed in previous studies based on a 
smaller path length (s ee |Prochter et al. 2006b). Our result and 
that of Procht eFet al.| ( |2006b| l are different at less than 2cr. The 
present result is statistically more significant. 

Possible explanations of this excess include: dust obscura- 
tion that could yield such lines of sight to be missed in quasar 
studies; difference of the beam sizes of the two types of back- 
ground sources; gravitational lensing. In order to retrieve more 
information to test these hypotheses we investigate in detail the 
properties of the strong Mg n systems observed with UVES. We 
find that the equivalent width distribution of Mg n systems is 
similar in GRBs and QSOs. This suggests that the absorbers are 
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Table 8. Properties of the foreground damped Ly-a systems detected along the UVES GRB spectra. 







z 

(Sub-)DLA 


(N in cm 2 ) 


log Noi log A'sui 
(N in cm" 2 ) (N in cm" 2 ) 


log Wpeii 

(Nincnr 2 ) 


[Fe/Si] [O/H] 


[Si/H] 


[Fe/H] 


GRB021004 


0.487 


There are no 


intervening (Sub-)DLA along this los 










(jKBOjO/JU 


0.542 


3.5655 


20.2 ±0.1 


N/A b < 14.7 


< 13.6 




< — 1.0 


< —2.3 


GRB 050s 20 A 


0.529 


2.3598 


20.1 ±0.2 


N/A b 13.84 ±0.02 


14.1 1 ± 0.02 


-0.22 ± 0.02 


-1,5 ± 0.2 


-1.7 ± 0.2 


GRB050922C 


0.468 


2.0778 


20.65 ±0.15 


> 15.30 14.38 ±0.10 


1 A A A i f\ f\*} 

14.44 ± 0.03 


0.11 ± 0.10 > -2.1 


1 Ol i A I O 

-1.83 ±0.18 


-1.72 ± 0.15 


GRB060607A 


0.596 


2.9374 c 


19.4 ±0.1 


> 15.07 > 14.74'' 


14.16 ± 0.03 


< -0.53 > -1.1 


> -0.2 


-0.7 ±0.1 






2.9322 c 


19.0 ±0.1 


< 12.8 < 12.6 


< 12.6 


< -2.9 


<-1.9 


< -1.9 


GRB071031 


0.540 


There are no intervening (Sub-)DLA along this los 










GRB080310 


0.501 


There are no intervening (Sub-)DLA along this los 










GRB080413A 


0.503 


There are no intervening (Sub-)DLA along this los 











a Redshift path calculated from the Ly-fl GRB absorption line to 3000 km/s from the GRB Ly-a, except for GRB 050730 where the z range ends at the beginning of the spectral gap at 
A = 5096A. The line is blended with other lines in the Ly-a forest. c To calculate the DLA number density we counted the z = 2.9374 and z = 2.9322 systems as a single system. 
d Blended with the Si iv/11402 absorption line associated to the GRB. 



more extended than the beam size of the sources which should 
not be the case for the different beam sizes to play a role in ex- 
plaining the excess (Por ciani et al.|2007] >. In addition, we divide 
our sample in three redshift bins and we find that the number 
density of strong Mg n systems is larger in the low redshift bins, 
favoring the idea that current sample of GRB lines of sight can be 
biased by gravitational lensing effect. We also estimate the dust 
extinction associated to the strong GRB Mg n systems and we 
find that it is consistent with what is observed in standard (sub)- 
DLAs. It therefore seems that the dust-bias explanation has little 
grounds. 

We tentatively infer the H i column densities of the strong 
systems and we note that the number density of DLAs per unit 
redshift in the UVES sample is probably twice larger than what 
is expected from QSO sightlines. As these sytems are expected 
to arise from the central part of massive haloes, this further sup- 
ports the idea of a gravitational lensing amplification bias. This 
hypothesis could be also supported by the results recently ob- 



tained by (Chen et al. 2008). These authors analyzed 7 GRB 
fields and found the presence of at least one addition object at 
angular separation from the GRB afterglow position in all the 
four fields of GRBs with known intervening strong Mgll galax- 
ies. In contrast, none is seen at this small angular separation in 
fields without known Mg II absorbers. 

We searched the Lyman-Q' forest probed by the UVES spec- 
tra for the presence of damped Lyman-ff absorption lines. We 
found four sub-DLAs with log N(H i) > 19.3 over a redshift 
range of z = 4.3. This is again twice larger than what is expected 
in QSOs. However the statistics is poor. It is intriguing that these 
systems are all located in the half redshift range of the Lyman-o- 
forest closest to the GRB. It is therefore not excluded that part 
of this gas is somehow associated with the GRBs. In that case, 
ejection velocities of the order of 10 to 25 000 km/s are required. 
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